A complete family of separability criteria 
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We introduce a new family of separability criteria that are based on the existence of extensions 
of a bipartite quantum state p to a larger number of parties satisfying certain symmetry properties. 
It can be easily shown that all separable states have the required extensions, so the non-existence 
of such an extension for a particular state implies that the state is entangled. One of the main 
advantages of this approach is that searching for the extension can be cast as a convex optimization 
problem known as a semidefinite program (SDP). Whenever an extension does not exist, the dual 
optimization constructs an explicit entanglement witness for the particular state. These separability 
tests can be ordered in a hierarchical structure whose first step corresponds to the well-known 
Positive Partial Transpose (Peres-Horodecki) criterion, and each test in the hierarchy is at least 
as powerful as the preceding one. This hierarchy is complete, in the sense that any entangled 
state is guaranteed to fail a test at some finite point in the hierarchy, thus showing it is entangled. 
The entanglement witnesses corresponding to each step of the hierarchy have well-defined and very 
interesting algebraic properties that in turn allow for a characterization of the interior of the set of 
positive maps. Coupled with some recent results on the computational complexity of the separability 
problem, which has been shown to be NP-hard, this hierarchy of tests gives a complete and also 
computationally and theoretically appealing characterization of mixed bipartite entangled states. 



I. INTRODUCTION 

Entanglement is one of the most fascinating features 
of quantum mechanics. As Einstein, Podolsky and 
Rosen Ij pointed out, the quantum states of two phys- 
ically separated systems that interacted in the past can 
defy our intuitions about the outcomes of local measure- 
ments. Entangled pure states have zero entropy but can 
appear to have maximal entropy when the experimenter 
only has access to one of the subsystems. On the other 
hand, Bell inequalities quantify the extent to which 
local measurements on separated quantum systems can 
be correlated in ways that are forbidden in any local clas- 
sical model. Violations of these inequalities require en- 
tanglement. Moreover, it has recently been recognized 
that entanglement is a very important resource in quan- 
tum information processing, allowing certain important 
tasks such as teleportation, quantum computation, quan- 
tum cryptography and quantum communication to name 
a few Q. 

For the case of pure states, determining when a given 
state is entangled is very easy, since it is based on prop- 
erties of the Schmidt decomposition or, equivalently, the 
rank of the reduced density matrices, which can be com- 
puted very efficiently. However, for the case of mixed bi- 
partite states, no single practical procedure that can be 
guaranteed to detect the entanglement of every entangled 
state has been found. In the past few years, a consider- 
able effort has been dedicated to this problem 0,0,IEI3- 
Still only incomplete criteria have been proposed, that 
can detect some entangled states but not all of them or 
that work only for certain restricted dimensions. This is 
a somewhat uncomfortable situation, since all the quan- 
tum states generated in the laboratory for practical ap- 
plications of quantum information processing are mixed 



states. Hence the need, not only from the theoretical but 
also from the practical point of view, of having an effi- 
cient tool that would allow us to determine when a given 
state is entangled. 

A bipartite mixed state is said to be separable |^ (not 
entangled) if it can be written as a convex combination 
of pure product states 



P = ^Pt\'>Pi){'^h\ ® \(t)t){4>t\ 



(1) 



where jV'i) and are state- vectors on the spaces Ha 
and TLb of subsystems A and B respectively, and pi > 
0, Pi = 1. If a state admits such a decomposition, then 
it can be created by local operations (unitary transforma- 
tions, measurements, etc.) and classical communication 
(LOCC) by the two parties, and hence it cannot be an 
entangled state. Despite the simplicity of (^3), it has been 
shown recently by Gurvits that deciding whether or 
not such a decomposition exists for a given density ma- 
trix is an NP-hard problem. This result destroys any 
hope of finding a computationally efhcient tool to de- 
termine entanglement of mixed states as was the case for 
pure states, so long as the widely believed result P ^ NP 
is actually true. But there are some instances of the sep- 
arability problem that allow efficient algorithms to solve 
them. This is one of the basic ideas behind separability 
criteria. 

A separability criterion is based on a simple property 
that can be shown to hold for every separable state. They 
provide necessary but not sufficient conditions for separa- 
bility. If some state p does not satisfy the property, then 
it must be entangled. But if it does satisfy it, that does 
not imply that the state is separable. One of the first and 
most widely used of these criteria is the Positive Partial 
Transpose (PPT) criterion, introduced by Peres 0. If 
p has matrix elements pik,ji = {i\ <8) {k\p\j) ® \l) then the 
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partial transpose p is defined by 

Pjkji = Pjk,ii ■ (2) 

If a state is separable, then it must have a positive partial 
transpose (PPT). To see this consider the decomposition 
(Q) for p. Partial transposition takes to | V'i' ) (V'l* I j 

so the partial tranpose of p can be written as 

Clearly p-^-* is a valid quantum state and in particu- 
lar it must be positive semidefinite. Thus any state for 
which p^^ is not positive semidefinite is necessarily en- 
tangled. This criterion is computationally very easy to 
check. Furthermore, it was shown by the Horodeckis 
based on previous work by Woronowicz 12], to be both 
necessary and sufficient for separability in 7^2 ® 'H2 and 
7^2 'Si 7^3. However, in higher dimensions, there are PPT 
states that are nonetheless entangled, as was first shown 
in , again based on ■ These states are called bound 
entangled states because they have the peculiar property 
that no ent ang lement can be distilled from them by local 
operations 

A different useful separability criterion, that has been 
used to show entanglement of PPT states, is the range 
criterion 0, ^| . It is based on the fact that for every 
separable state p there exist a set of pure product states 
{|'0i)|</'i)} that span the range of p while {\'fpi)\(t>i)} span 
the range of p^^ , as can be easily seen by looking at equa- 
tions and ||2J). This criterion is sometimes stronger 
than PPT, but in some cases it can be weaker (for exam- 
ple, when considering full rank non-PPT states).. Other 
criteria, that are in general weaker than PPT are the 
reduction criterion [TM IT^ and the majorization crite- 
rion ,17i] . None of these criteria, nor a combination of 
them are sufficient to give a complete characterization of 
separable states. 

Another approach to distinguishing separable and en- 
tangled states involves the so called entanglement wit- 
nesses (EW) 0. An EW is an observable W whose 
expectation value is nonnegative on any separable state, 
but strictly negative on an entangled state p. We say in 
this case that W "witnesses" the entanglement of p. Be- 
sides giving another theoretical tool to detect entangled 
states, this idea addresses the question of whether there is 
an experimental way of distinguishing an entangled state 
from a separable one. By studying the geometrical struc- 
ture of the set of quantum states, it can be shown that 
for every entangled state there exists an entanglement 
witness W 0,^2- Thus, there is always an observable 
that can be measured that will show that the state is 
entangled. 

There are two other important mathematical objects 
related to entanglement witnesses. Although these do 
not have the physical interpretation of observables they 
allow connections to other results in the mathematical 
and mathematical physics literature. In the first place. 



there is a correspondence that relates entanglement wit- 
nesses to linear positive (but not completely positive) 
maps from operators on Ha to operators on Hb (or vice 
versa); see equation (|89|l and reference Applying 
such a map to one half of an entangled state does not 
necessarily result in a positive matrix. For this reason 
positive maps were rejected as possible physical evolu- 
tions of quantum states in favor of the completely pos- 
itive maps. The PPT test has this structure where the 
transpose is the positive map. Any positive but not com- 
pletely positive map results in an analogous separability 
criterion. The equivalence between entanglement wit- 
nesses and positive maps implies that if p is entangled 
there is always a positive map that will detect the en- 
tanglement in this way [T^, [Q ■ The characterization of 
positive linear maps was in fact the original motivation 
for studying the separability question jl2l |. 

Finally, there is a well known mapping between pos- 
itive linear maps and positive semidefinite biquadratic 
forms |2oL I2H . This can be appreciated simply by writ- 
ing the condition that W is positive on pure product 
states explicitly in terms of the elements of W and the 
state vectors for the two systems, as in equation (|43ll in 
Section IVII This suggests the use of results from real 
algebraic geometry (see for example |2^ and the refer- 
ences therein) to attack the separability problem. Indeed, 
the semidefinite programming techniques we employ here 
were first developed in this general context [23l |. 

The question of whether a given state p is separable 
may be phrased as quantified polynomial inequalities in 
a finite number of variables: 

VW^[V|?A)V|(/)) {^\{(l)\W\^p)\(j)) > =^ Tt[pW] > 0]. (4) 

If this proposition is satisfied then p is separable. Since 
the inequalities may be expressed in terms of polynomi- 
als of the variables (the components of W, \ , this 
is a semi-algebraic problem. Much is known about the 
general class of semi- algebraic problems, in particular the 
fact that they are decidable. The Tarski-Seidenberg de- 
cision procedure [23 | can then be used to provide an ex- 
plicit algorithm to solve the separability problem in all 
cases and therefore to decide whether p is entangled. A 
drawback of this approach is that most exact techniques 
in algebraic geometry scale very poorly with the number 
of variables (the Hilbert space dimensions in the separa- 
bility problem). For this problem, these methods do not 
perform well in practice except for very small problem 
instances. This is in contrast to the PPT test which may 
be implemented very efficiently but does not always settle 
the question of separability of p. In this paper we discuss 
a set of separability criteria that also have this prop- 
erty; they all scale polynomially with the Hilbert space 
dimension and perform well in practice, any state p that 
is entangled is detected by one of the tests but no one 
test detects all entangled states. Since the separability 
problem is NP-Hard it is very unlikely that a procedure 
guaranteed to solve the problem in all instances can scale 
well with Hilbert space dimension. As a result our fam- 
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ily of separability tests is, in some sense, the best way 
of solving the problem from a practical point of view, in 
that simple tests will detect the easiest instances of the 
problem, while the more complicated instances genuinely 
require more computational resources. 

The most important characteristic of the separability 
problem is the fact that the separable states form a con- 
vex set. The existence of entanglement witnesses, ob- 
servables that are positive on separable states but neg- 
ative on some entangled state, is a direct result of this 
convexity. There has been much work on the separabil- 
ity problem, particularly from the Innsbruck-Hannover 
group as reviewed in 0, 0|, that emphasizes convexity 
and proceeds by characterizing entanglement witnesses 
in terms of their extreme points, the so-called optimal 
entanglement witnesses, and PPT entangled states in 
terms of their extreme points, the edge PPT entangled 
states 0, '25^ . Convexity also plays a central role in our 
work which provides a computational means of construct- 
ing entanglement witnesses with certain properties. It is 
interesting that our construction will allow us to charac- 
terize the interior of the set of entanglement witnesses, 
but not its extreme points. 

Beyond the separability problem, many problems of 
interest in quantum information have the structure of 
convex optimizations j^, a fact that has found increas- 
ing application in the field in recent years. One early 
example is the use of results about linear programming 
to find the optimal local entanglement concentration pro- 
cedure for a pure bipartite state in Our work will 
involve convex optimizations known as semidefinite pro- 
grams [2^ generalizations of linear programs that 
optimize a linear function of a positive matrix subject 
to linear constraints. Semidefinite programming argu- 
ments have also been used been used in the quantum in- 
formation literature to address questions about quantum 
coin tossing, distillation and optimal state transforma- 
tions 

In this paper we discuss in detail a family of separabil- 
ity criteria introduced in [s^, that can be ordered into 
a hierarchy of tests that have the following two very im- 
portant properties: i) the hierarchy is complete, i.e., any 
entangled state will be detected by some test in the hi- 
erarchy, ii) there are efficient computational algorithms 
to check each of the tests. This provides us with a very 
practical algorithmic way for testing entanglement of a 
given bipartite mixed states, that is guaranteed to detect 
any entangled state. Furthermore, the algorithm con- 
structs an explicit proof of this fact in the form of an 
entanglement witness. This in turn helps us to develop a 
characterization of almost all positive maps that are not 
completely positive. 

The paper is organized as follows. In Section II we in- 
troduce a new family of separability criteria; in Section 
III we introduce and discuss the properties of semidefinite 
programs (SDP), and show that each separability test in 
the family can be cast as a SDP, and briefly discuss the 
resources needed to implement them. In Section IV we 



discuss how to take advantage of the symmetries that 
each test requires to further reduce the computational 
resources needed. In Section V we present an explicit 
proof of the completeness of the hierarchy, translating 
previous results 35. 36J into the language of density ma- 
trices. Section VI shows how the duality of the SDP can 
be exploited to construct an entanglement witness that 
proves entanglement for a given state, and we discuss the 
algebraic properties of these witnesses. In Section VII we 
present examples of the application of the hierarchy. In 
section VIII we discuss how to construct an entangled 
state that is not detected by the second test of the hi- 
erarchy and present several important consequences of 
this result. Section IX shows how to use an SDP to test 
indecomposability of an entanglement witness and con- 
struct a bound entangled state detected by it. In Section 
X we discuss the connection betweeen entanglement wit- 
nesses and positive maps, and show how the properties of 
the witnesses obtained through our hierarchy of tests can 
be translated into a characterization of strictly positive 
maps. Finally, in Section XI, we summarize our results 
and present our conclusions. 

II. PPT SYMMETRIC EXTENSIONS AND 
NEW SEPARABILITY CRITERIA 

Any separable state p in Ti.A (E)Ti.B can be written as in 
Consider now the state p in Ha ^H-b ® Ha, given 

by 

p = \<l>^){<l>^\ ® (5) 

Then p has the following properties: i) p is an extension 
of p to three parties, in the sense that 

Trc[p]-P, (6) 

where Trc means that we take the partial trace over the 
third party which we have taken to be equal to Ha- h) 
p is symmetric under interchanges of the first and third 
parties, i.e., the two copies of party A. More precisely, if 
we define the swap operator P by 

P\t)®\k)®\j) = \j)®\k)®\i), (7) 

the symmetry condition can be written as 

p = PpP. (8) 

iii) p must remain positive under any partial transposi- 
tion (since p is also a separable state). Note that, due 
to the symmetry taking partial transpose with re- 
spect to the third subsystem is equal to taking it with 
respect to the first one. Now, for an arbitrary state p in 
Ha^Hb, we will call p a PPT symmetric extension of p 
to two copies of Ha , if and only if, p satisfies the three 
properties stated above. Since we have shown by con- 
struction that any separable state has a PPT symmetric 
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extension to two copies of Ha, then we can use its exis- 
tence as a separabihty criterion. If a given state does not 
have such an extension, then the state must necessarily 
be entangled. 

We can take this idea of the existence of PPT symmet- 
ric extensions further by considering extending the state 
to an arbitrary number of copies of subsystem A. For any 
separable state in T-La (^Hb given by ^ , the state 

p = ® ® (9) 

is a state in 1-La®'Hb®'H^^'^ that, i) is symmetric under 
interchanges of any two copies of subsystem A, ii) yields 
the original state p in 'Ha®'Hb when we trace out any 71— 
1 copies of subsystem A, and iii) remains positive under 
all possible partial transpositions. Again, for an arbitrary 
state p, we will call p a PPT symmetric extension of p to 
n copies of party A, if and only if, p satisfies properties i), 
ii) and iii). And as before, we can use the existence of this 
extension to n copies of subsystem A as a separability 
criterion. We have thus generated a countably infinite 
family of separability criteria. Note that the same idea 
can be generalized to the multipartite case: the existence 
of PPT symmetric extensions to any number of copies of 
the parties is a separability criterion. 

For the bipartite case, these separability criteria are 
not completely independent of each other, but they ac- 
tually have a hierarchical structure. We will now show 
that if a state has a PPT symmetric extension to n copies 
of A, call it pn , then it must have a PPT symmetric ex- 
tension to n — 1 copies of A. Let pn-i = Tr^[jO„], where 
A represents one of the copies of A. It is easy to see 
that Pn-i will inherit from p„ the property of being sym- 
metric under interchanges of copies of party A, since we 
have just removed one of the copies. It is also obvious 
that Pn-i is an extension of p to n — 1 copies of A. Let's 
assume that is not PPT. Then there is a subset T of the 
parties such that p^L^ ^ negative eigenvalue, where 
Ti represents the partial transpose with respect to all the 
parties in subset J. Let |e) be the corresponding eigen- 
vector and let be a basis of the system A over which 
the partial trace was performed. Since p„ is PPT, then 
{e\{i\p^^\e)\i) > 0, for aU i. Then 

Y.{e\{^\pl^\e)\^) = {e\TTA[pl^]\e) >0. (10) 

i 

Since we performed the partial trace over a party that is 
not included in X , we can commute the trace and the 
partial transpose, and using pn-i = Tr aIpu], we have 
{e\Pnti\e) > 0, which contradicts the fact that |e) is an 
eigenvector of pJ{Li with negative eigenvalue. 

We have then constructed a family of separability cri- 
teria with a natural hierarchical structure. If we take the 
usual PPT criterion as the first step of the hierarchy, the 
existence of a PPT extension to two copies of A as the 
second step, and so on, we see that the tests are ordered 
in such a way that each test is at least as powerful as the 



previous one, in the sense that if a state was shown to 
be entangled by one of them, it will be also shown to be 
entangled by all the tests that are higher in the hierar- 
chy. This family of tests has several very important and 
useful properties. It can be shown that each test can be 
cast as a semidefinite program (SDP), which is a class 
of convex optimization problems for which efhcient algo- 
rithms exist. The duality structure of the SDP allows us 
to construct an explicit entanglement witness whenever 
a state fails one of the tests. And finally, it can be proven 
that the hierarchy is complete, i.e., every entangled state 
is guaranteed to fail the test at some finite point in the 
hierarchy. 

III. SEMIDEFINITE PROGRAMS AND 
SEARCHING FOR PPT SYMMETRIC 
EXTENSIONS 

In this section we will introduce and discuss the struc- 
ture of a semidefinite program and we will show explicitly 
how to apply it to the problem of searching for a PPT 
symmetric extension. 

A. Semidefinite programs 

A semidefinite program (SDP) is a particular type of 
convex optimization problem [29 . l28l |. An SDP corre- 
sponds to the optimization of a linear function subject 
to a linear matrix inequality (LMI). A typical SDP has 
the form 

minimize c^x 

subject to F(x) > 0, (11) 

where c is a given vector, x = (a;i, . . . , a;„), and -F(x) = 
_Fo -t- XiFi, for some fixed hermitian matrices Fi. The 
inequality in the second line means that the matrix -F(x) 
must be positive semidefinite. The minimization is per- 
formed over the vector x, whose components are the vari- 
ables of the problem. The set of feasible solutions, i.e., 
the set of x that satisfy the LMI, is a convex set. In the 
particular case in which c = 0, there is no function to 
minimize and the problem reduces to whether or not the 
LMI can be satisfied for some value of the vector x. In 
this case, the SDP is refered to as a feasibility problem. 
The convexity of the SDP has made it possible to de- 
velop sophisticated and reliable analytical and numerical 
methods to solve them '28]. 

A very important property of a SDP, both from the 
theoretical and applied points of view, is its duality struc- 
ture. To any SDP of the form (|ll|l . which is usually 
called the primal problem, there is associated another 
SDP, called the dual problem, that can be stated as 

maximize — Tr[fo-^] 
subject to Z > 

Tt[F,Z] = c„ (12) 



5 



where the matrix Z is hermitian and is the variable over 
which the maximization is performed. This corresponds 
to the maximization of a hnear functional, subject to 
linear constraints and a LMI. Let x and Z be any two 
feasible solutions of the primal and dual problems respec- 
tively. Then we have the following relationship 

c^x + Tr[FoZ] = Tr[F(x)Z] > 0, (13) 

where the last inequality follows from from the fact that 
both i^(x) and Z are positive semidcfinitc. From and 
(|12|l we can see that the left hand side of (|13|l is just the 
difference between to objective functions of the primal 
and dual problem. The inequality in H13|) tells us that the 
value of the primal objective function evaluated on any 
feasible vector x, is always greater or equal than the value 
of the dual objective function evaluated on any feasible 
matrix Z . This property is known as weak duality. Thus, 
we can use any feasible x to compute an upper bound 
for the optimum of —Tt[FoZ], and we can also use any 
feasible Z to compute a lower bound for the optimum of 

T 
C X. 

If the feasibility constraints on both the primal and 
dual problems are satisfied for some Z > and x such 
that -F'(x) > 0, the problems are termed strictly feasi- 
ble, and the optimum values of the primal and dual for- 
mulations are equal. This property is called strong du- 
ality. Furthermore, there is a feasible pair (xopt,-^opt) 
achieving the optimum. In this case, as can be seen from 
equation (fT^ . we have Tr [F(xopt)^opt] = 0, and thus 
F(xopt).^opt = 0, so the hermitian matrices F(xopt) and 
Zopt have orthogonal ranges. This is known as the com- 
plementary slackness condition '28l . 

Equation H13|l has another important application. 
Consider the particular case of a feasibility problem (i.e., 
c = 0). Then, equation H13|l will read 

Tr[FoZ] > 0, (14) 

and this must hold for any feasible solution of the dual 
problem. This property can be used to give a certificate of 
infeasibility for the primal problem: if there exists Z such 
that Z > and TT[FiZ] = 0, that satisfies Tr[_Fo-^] < 
0, then the primal problem must be infeasible. We will 
show later that for the particular case of our hierarchy of 
separability tests, whenever a PPT symmetric extension 
of p cannot be found (primal problem is infeasible), the 
certificate provided by the dual problem is nothing but 
an entanglement witness for the state p. 

B. Separability tests as semidefinite programs 

Each test in the hierarchy of separability criteria in- 
troduced in Section II can be written as a semidefinite 
program. We will show in detail how the SDP is setup 
for the second test in the hierarchy, which corresponds 
to searching for PPT symmetric extensions of p to two 
copies of subsystem A. The general case, of extensions to 
n copies of party A, can be constructed in a similar way. 



Let {crfj.^i, {cTf}j=i be bases for the space of hermi- 
tian matrices that operate on and Ti^ , of dimensions 
dA and ds respectively, such that they satisfy 

Ti[afaf]^axS,j and Tr[(Tf ] = (5,i, (15) 

where X stands for Aoi B, and ax is some constant — the 
generators of SU(n) could be used to form such a basis. 
Then we can expand p in the basis {af ® <^f}^ and write 

P = Ey P^3'^f «> '^f , with = a^'ag^Tvipaf (g) af ]. 
In the same way, we can expand the extension p in Ha <8) 
Hb (8) Ti-A as 

P - "^^fc ^'^^^} + 

i<k 
kj 

where we made explicit use of the swapping symmetry 
between the first and third parties, that we require from 
p. To satisfy the condition that p is an extension of p, 
we need to impose 

Trc[p] = p, (17) 

where Trc means tracing out the third party. Using 115() 
and the fact that {af" ® a^} form a basis of Ha ® Hb, 
equation l(T7|) reduces to 

Piji = Pij- (18) 

This fixes some of the components of p. The remaining 
ones will play the role of the variables in the SDP. The 
LMIs come from requiring that the extension p and all its 
partial tranposes are positive semidefinite. If we define 

Go = X! ® '^f + 

j 

G^Ji = af ®af ®af i>2, 

Gijk = [crt ® erf ®^k +(^k® erf ®af) k> i> 2, 

(19) 

we can write the PSD condition p > as 

G(x) = Go + ^ xjGj > 0, (20) 
J 

where we have collected all the subindices in H19|l into 
one subindex J. Equation H2UI) has exactly the form that 
appears in (|ll|l . The role of the variable x is played by 
the coefficients pijk{k ^ l,k > i), which can vary freely 
without affecting the extension condition (|17|l . The num- 
ber of free variables is m = {d\d^ — d'^d^)/2, and each 
matrix Gj has dimension n = d\dB. Since p is symmet- 
ric under swaps of the first and third parties, there are 
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only two independent partial transpositions that can be 
applied to it, which we can take as partial transposes with 
respect to the first and second parties (one of the copies 
of A, and subsystem B). The requirement that these two 
partial transposes are positive leads to two more LMIs, 
given by 



P^-* > 



and 



>o, 



(21) 



where the Gj matrices for these two inequalities are re- 
lated to the matrices given in l|19() by the appropriate 
partial transposes, namely GJ-* and G^^. We can ac- 
tually combine the three LMIs into one, by defining the 
matrix 



F 



(22) 



and using the fact that a block-diagonal matrix C = 
A(B B is positive semidefinite, if and only if, both A and 
B are positive semidefinite. 

We have then stated the search of a PPT symmetric 
extension of p as an SDP, in which the objective function 
is zero (c = 0), so it corresponds to a feasibility problem, 
and the LMI condition reads F = p® p^^ © p^^ > 0, 
which encodes the requirement of the extension and its 
partial transposes being positive semidefinite. The SDP 
will then take the form 



minimize 
subject to F >Q. 



(23) 



In appendix A, we discuss a slightly modified version of 
the SDP that has the advantage of performing better 
numerically, but we will keep the form (|23() for all the 
analytical discussions, since its dual form is more clearly 
related to the construction of entanglement witnesses, 
which is another one of the main results of this paper. 

The SDP for other tests of the hierarchy (extensions to 
more copies of party A), can be constructed in the same 
way, by generating the matrices Gj with the appropriate 
symmetry, and constructing the block-diagonal matrix _F, 
whose blocks correspond to all the independent partial 
transposes that can be applied to the extension p. 



C. Resources needed to implement the tests 

As we mentioned before, there are very efficient algo- 
rithms to solve semidefinite programs, and we can use 
their properties to discuss the computational resources 
required to implement a general step in our hierarchy 
of tests. Assume that we are searching for a PPT sym- 
metric extension of a state p in T-La (8) Tis to fc copies of 
subsystem A, with cIa and the dimensions of T-La and 
T-Lb respectively. Then, the corresponding semidefinite 



d\ variables 



program will have m 

and a matrix G with (fc -f 1) blocks of dimension (P^dn 
Numerical SDP solvers are described in detail in [23 



Typically they involve the solution of a series of least 
squares problems each requiring a number of operations 
scaling with problem size as 0{m^n^), where F{x) is an 
nxn matrix. For SDPs with a block structure these break 
into independent parts each with a value of n determined 
by the block size. The number of iterations required is 
known to scale no worse than 0{n^^^). Thus for any fixed 
value of k the computation involved in checking our crite- 
ria scales no worse than 0{(f^''^^) which is polynomial in 
the system size. On the other hand, for d-A and ds fixed, 
the size of the matrix F(x) scales exponentially with the 
number of copies k. There is, however, a significant im- 
provement that can be accomplished by exploiting the 
swapping symmetry to its fullest. In the next section we 
will show that we can impose a stronger restriction on 
the extension that brings the scaling of resources down 
to polynomial in the number of copies of subsystem A, 
for fixed dA, ds- It is important to point out that the 
resources required to solve the separability problem have 
been proven to scale super polynomially only when the 
dimensions of both systems are allowed to vary. These 
two results are consistent, although the complexity re- 
sult implies that there is no value of k such that the A:*'* 
test detects all entangled states for all values of dA, ds- 



IV. EXPLOITING THE SYMMETRY 

Each test for separability searches for an extension of 
a state p in TLa ® T^b of dimension dAds, to the space 
nf'^'HB, that has dimension d\dB (where, without loss 
of generality, we have interchanged the order of Tis and 
all copies of A for convenience). We see that the dimen- 
sion of the extended space increases exponentially with 
the number of copies of party A. We have shown that we 
can impose further restrictions on the extension, and in 
particular we require it to be invariant under swaps of the 
copies of subsystem A. This reduces the size of the space 
over which we search for the extension, but the scaling 
with the number of copies remains exponential, which is 
not desirable of a practical tool for deciding separability 
of a state. However, we can actually impose a stronger 
constraint on the form of the extension, that reduces the 
scaling of its size from exponential to polynomial in the 
number of copies. 

As we pointed out before, any separable state in Ha ^ 
He of the form Q has a PPT symmetric extension to 
Tif^ i^TCb, that we can explicitly write as 



(24) 



This extension is obviously invariant under swaps of 
copies of A, and we used this property to restrict the 
form of the matrices Fj in the LMI of our SDP. But p 
has a more constraining property: its support and range 
are contained in the symmetric subspace of ® TCb 
(where the symmetry is understood to apply only to the 
copies of A). For the case of the extension to two copies 
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of system A, we can write the projector into this symmet- 
ric subspace as tt = ^(1 + P), with P the swap operator 
defined in ((JJ). Then, the symmetry requirement on the 
extension takes the form p = npi:. 

For an arbitrary p, we can now restrict our search to 
extensions that satisfy this property. If {Sf} is a basis 
of hermitian matrices having support and range in the 
symmetric subspace of Ti®'^, this restriction is equivalent 
to only considering matrices G in 1)201) of the form G — 
sf (8) . Since the dimension of the symmetric subspace 

in nf" is 



ds. = \' - , (25) 



with dA the dimension of TY^i, the number of matrices of 
this form is dig^d"^. The number of variables in the SDP 
is this number minus the number of constraints given by 
equation H17() . which is m = {d'g^ —d\)d^. By using (|25|l . 
we get m ~ 

0{k'^dA-i)^^ which for a fixed size of A is only 
polynomial in the number of copies. Since the matrices 
Gj have range and support only on this symmetric sub- 
space, we know that by a suitable change of basis they 
can be simultaneously block-diagonalized, with the only 
nonzero block having size n — d^^d^. 

The SDP that searches for the PPT extension also re- 
quires to check positivity of a certain number of partial 
transposes of the extension p. These checks translate into 
a bigger LMI, although we will now show that this does 
not change the scaling properties of its size. Consider 
the case in which we apply the partial transpose to the 
first I copies of A, which we will denote p^A»' . Since the 
matrices Gj have support and range only on the sym- 
metric subspace of Ti^*^, it is not difficult to show that 

the matrices Gj"*®' must have support only on the ten- 
sor product of a subspace isomorphic to the symmetric 
subspace o/7i^' and the symmetric subspace ofTi^^ 
The dimension of this tensor product is just the product 
of the dimensions of the two subspaces. If we perform a 
change of coordinates by rotating to a basis that contains 
a basis of this tensor product of symetric subspaces, we 
can see that the size of the matrices G^^®' can be taken 
as d'g^d'g^^ This scales at most as 0{k'^^''"^~^^). Since 
the number of independent partial tranposes is (fc -f- 1), 
as a result of the symmetry requirements, the size n of 
the matrices in the LMI scales not worse than 0{k^'^^^^). 
Combining this with the scaling of the number of vari- 
ables m shown above and the scaling properties of solving 
the SDP, which is given by O(m^n^), we can see that for 
fixed dA, the tests in the hierarchy scale as 0(fc(^'' -*""')), 
which is polynomial on the number of copies of party A 
for a fixed dA- 



V. COMPLETENESS OF THE HIERARCHY OF 
TESTS 

One of the main results of this paper is the complete- 
ness of the hierarchy of separability tests. This result 
allows us to give an algorithm that will show if a state 
is entangled in a finite number of steps (although this 
number may be high for some states). Even though the 
hierarchy of tests is a new result, the proof of its com- 
pleteness is identical to the proof of certain properties of 
the possible equilibrium states of a system that interacts 
with a thermal bath. These results, which were proved 
by Raggio et al. [s^ , and Pannes et al. [s^ , have been in 
the literature for quite some time. 

It was noted in [s^l that this result |3^ [s^ could be 
interpreted as a characterization of bipartite quantum 
states, that requires that the only states that can have 
symmetric extensions to any number of copies of one of its 
subsystyems, are the separable states. The same idea was 
independently conjectured recently by Schumacher [s^ . 

We will present a proof of the completeness of the hier- 
archy, which is basically the proof found in _35j, applied 
to the case of bipartite mixed states on finite dimensional 
spaces. Our discussion has the same level of mathemati- 
cal rigor and is based on the techniques presented in the 
discussion of the Quantum de Finetti Theorem in |39j |. 
The theorem we will prove is stronger than our hierarchy, 
since it requires the existence of symmetric extensions 
without any requirements on the partial transposes. The 
completeness of our hierarchy can be deduced from this 
result as a corollary. 

Theorem 1 (Completeness) Let p be a bipartite 
mixed state in Ha Ti-B ■ Then p has a symmetric exten- 
sion to k copies of subsystem A for any k, if and only if, 
p is separable. 

Proof: One of the implications is trivial. Assume p is 
separable. Then we can write 

p = (26) 

^From this expression, we can write down explicitly a 
symmetric extension p for any value of fc, namely 

p = (27) 

and this completes the first part of the proof. 

To prove the other implication, the idea is to use the 
existence of the extensions to construct a set of states 
in Ha"' that can be shown to be separable by using the 
Quantum de Finetti Theorem, and then show that this 
result implies that the extensions themselves have to be 
separable. Let p be a state in Ha ^ Hb such that for 
any n, there is a symmetric extension of p in 7i®" ® Hb, 
which we will call /5„. Let us pick a fixed value k for 

the number of copies of party A. Let the set be 
a basis for the set of hermitian operators in Hb, such 
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that bi > for all i (i.e., all these operators are positive 
definite j40]), and in particular let us choose bi = 1^, 
the identity in Hb- Now we define the operator 



Pbi,k 



<^ bi)pk\ 



(28) 



where 1^ is the identity on subsystem A. The operator 
Pbi,k is positive semidefinite (PSD) and non zero since 
all the operators bi were taken to be strictly positive. 
Then pt^^k is proportional to a state in TiJ^ , since it is 
hermitian and PSD. We can choose the operators bi such 
that Tr[pb;jt] = 1 for all fc, so that (|28|l is actually a 
normalized state in Ti.'^ . 

We will now prove that the existence of symmetric ex- 
tensions pk of p for all fc, imply that we can choose the 
states Pbi,k to be exchangeable 39]. Recalling the defi- 
nition of exchangeability we need to show that, for any 
I > 0, there are states Pbi,{k+i) that are symmetric and 
satisfy 



Pb„k = Tr^fc+i...Afe+,[p6i,(fc+i)]- 



(29) 



Let us fix k and assume that there is not an extension jjk 
such that the state pb^^k given by is exchangeable. 
That means that there has to be a value h for which 
equation (|29|l is not satisfied for any pb-^k and pbi,(fc+ii)- 
But since p has symmetric extensions for all k, we can 
just choose an extension to {k + li) copies p(^k+h)^ and we 
have that Tr^^^^...^^^,^ [P{k+h)] is a symmetric extension 
of p to A: copies, and 

Pb„k = Ttb [{If 'E>b,)TrA,^,...A,^,Jp(k+i,)]] 



= Tta 



k + l-'-^k + li 



TTB[{tT^'^^b,)p<^k+l.^] 



1 [Pb. 



{k+h)\ 



(30) 



This is a contradiction, so we can conclude that we can 
always choose the states pb^^k to be exchangeable. 

The state pb^^k satisfies then the hypothesis of the 
Quantum de Finetti Theorem [39| , and so we know there 
is a unique probability measure function Pb- (g) > 0, such 
that 



Pbi,k 



D 



g^^PbMdQ, 



(31) 



where D represents the space of states in Ha (i-e., the 
set of hermitian, positive semidefinite operators of trace 
one). 

For each g, we can think of Pbi{g) as a functional ap- 
plied to the operators 6^, which we will note Fg, defined 
as Fg(bi) = Pb-{g). This functional is linear on con- 
vex combinations of positive operators. To see this, let 
^ > 0. Then Fg{pbi + {1- p)bj) = Pf,b,+(i-ii)bj{Q), where 
^A'fci+(i-/i)&j is the unique probability density that satis- 



D 



fies 

P(/i6i + (l-Ai)6j),fc / Q^''Ptibi + {l-fi)b.iiQ)dg 

Jd 

= Tib [{tf (g) {pb, + (1 - p)b,))pk] 
= pTiB [{tf ® b,)pk] + 
+ {1-p)Ttb [{if ^b,)pk] 

{^,P,^{g) + {l-f,)P,^{g)) g^f^dg. 

(32) 

The second equality in (|32|l holds because we are con- 
sidering a convex combination of the operators bi , which 
guarantees that TrB[(l®''' (g) {pbi -|- (1 — p)bj))pk] is nor- 
malized. Then, by the uniqueness of the probability den- 
sity in the Quantum de Finetti Theorem, we have 

f'Mfci+(i-M)&j = pPbAo) + (1 - p)Pb,{g), (33) 
which translates into 

Fgipb^ + (f - P)bj) = pFg{bi) + (1 - p)Fg{b,). (34) 

Then Fg is a linear functional on convex combinations of 
positive states vhTLb- 

Since Fg is defined on a basis, there is a unique way of 
extending this functional linearly to the whole space of 
operators in Hb- So we have a linear, positive and con- 
tinuous functional on a finite dimensional Hilbert space, 
and it is a well-known result that any such functional can 
be written as 



Fg{b)^TvB[agb] V6, 



(35) 



for some unique positive semidefinite operator dg in TLb- 
This operator might not be a state in TLb since it need 
not be normalized. We can then define a function 



P{g) - Tr[a,], 



(36) 



that is nonnegative. If P{g) is nonzero, we can define 
ag — ag/P{g). Then H35|l takes the form 



Pb{g)=Fg{b)^TrB[(Jgb]P{g) V6. 



(37) 



Note that since dg is normalized, P{g) — P^^{g), which 
shows that P{g) is a probability density. Using (|37|) in 
(EU, we get 



Pb,.k 



g''^''TTB[<7MPi9)dg 



D 



Trf 



(if^b,) / g'"'(^agPig)dg 
Jd 



(38) 



If P(p) = for some g, we can define ag arbitrarily, since 
it would not contribute to the integral in . Since 
is valid for all the elements bi of a basis of hermitian ma- 
trices in Hb, by comparing the expression in the second 
line with (|28|l . we can deduce that 



Pk 



g®'' ®UgP{g)dg. 



D 



(39) 
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This means that pk is a separable state, since (I39|l is an 
explicit decomposition as a convex combination of prod- 
uct states. Furthermore, since pk is an extension of our 
original state p, we have 

p = TYA2...AdPk\ 

= / Q®(JgP{Q)dQ, (40) 
J D 

which shows that p has to be a separable state. This 
concludes the proof of the theorem. □ 

It is clear that this theorem implies the completeness 
of the hierarchy of separability tests introduced in Sec- 
tion II, since a state that has PPT symmetric extensions 
to k copies of party A for all values of k obviously has 
symmetric extensions for all values ofk, which according 
to the theorem implies that the state must be separable. 
However, it is interesting to note that the PPT require- 
ment is not essential for the completeness of the hier- 
archy. Searching just for symmetric extensions is also 
a complete family of separability criteria and one that 
requires less resources. 

In local hidden variable (LHV) theories were also 
constructed for quantum states possessing so-called sym- 
metric quasi-extensions, where rather than requiring that 
the extension be positive as a matrix it is only required 
that it is positive on product states. The number of ex- 
tensions corresponds to the number of independent local 
measurement settings that the theory is able to describe. 
In fact our argument that only separable states have an 
arbitrary number of symmetric extensions generalizes to 
this case. Essentially all that is needed is a version of 
the Quantum de Finetti theorem that holds for entangle- 
ment witnesses as well as states but it is straightforward 
to check that the argument of 39] holds in this case also 
since only POVMs that act as tensor products on each 
subsystem are used in the proof. Hence although the use 
of quasi-extensions is strictly stronger for a small num- 
ber of local measurement settings, if the LHV is required 
to work for an arbitrary number of local measurement 
settings, the construction will only work for separable 
states. 

We could generate more families of criteria, by search- 
ing for symmetric extensions that have to satisfy some 
other constraint, and this family of tests would still be 
complete although it would in general require more re- 
sources. If the constraint can be written in terms of 
linear equalities and LMIs, we could still use an SDP 
to implement the tests. Choosing between these many 
possibilities is a matter of how well they perform in ac- 
tual examples. It becomes a trade-off between how much 
more powerful the tests become when more constraints 
are placed on the extensions, and how much this increases 
the resources needed. Including the PPT requirement on 
the extension has the advantage that it guarantees that 
the second and higher tests in the hierarchy are stronger 
than the PPT criterion, and we have found this to be a 
good trade-off in practice. 



VI. CONSTRUCTION OF ENTANGLEMENT 
WITNESSES 

An entanglement witness (EW) for a state p is a her- 
mitian operator W that satisfies 

Tt[pW] < and Tr[p,epW^] > 0, (41) 

where psep is any separable state 0, 0|. It is clear 
that if H41|l is satisfied, then p cannot be separable, and 
W gives a proof of that fact. This property has a very 
nice geometric interpretation. Since the set of separable 
states is convex, any point that does not belong to it (like 
any entangled state) , can be separated from the set by an 
hyperplane. In our case, the operator W defines the hy- 
perplane. This result is known as the Hahn-Banach theo- 
rem ^]. In practice, finding a W satisfying Tr[pW^] < 
is not difficult, but proving Tr[psepW^] ^ might be very 
hard. To understand the reason for this, let us recall that 
any separable state can be written as a convex combina- 
tion of projectors into pure product states 

Psep 

= ^p,|a;)(x|®|2;)(2;|, (42) 

i 

where \x) = Yl,iXi\i), \y) = J^jUjU)^ for some bases 
{\i)} and {\j)} of Ha and Hb, respectively. Then, 
Tr[/OsepI^] > for any separable state Psep, if and 
only if, Tr[|a:)(a::| (g) > for any product state 

\x){x\ (8) \y){y\- Then we have 

Ewix,y)^ {xy\W\xy) = Tt[\x) {x\ ® \y) {y\W] 

= ^WijMX*y*Xkyi. (43) 

ijkl 

We can the interpret then the requirement that 
Tr [psepiy] > as a positivity condition on the biher- 
mitian form Ew{x,y) associated with the entanglement 
witness W ^ where bihermitian means that the form is 
hermitian with respect to x and hermitian with respect 
to y. It is a well-known result that checking positivity 
of an arbitrary real form is an NP-hard problem, and 
the result in !9| implies the same is true for bihermitian 
forms. This is the reason why constructing entanglement 
witnesses is not easy in general. 

As we mentioned in Section HI, any primal SDP has an 
associated dual problem that is also a SDP, and in partic- 
ular, whenever the primal problem is infeasible, the dual 
problem provides a certificate of this infeasibility. We 
will show that in the case of our separability tests, this 
infeasibility certificate generated by the dual problem is 
actually an entanglement witness. 

Consider the SDP and let us focus on the second 
test of the hierarchy, i.e., searching for PPT symmetric 
extensions to two copies of party A. In this case, the dual 
problem takes the form 

maximize — Tr[FoZ] 
subject to Z >0, 

Tt[E,Z] = 0, (44) 
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where Fq has three blocks that encode the extension and 
its two independent partial transposes, and from (|22|l we 
can see that it has the form 

Fo^Go®G^^ OG^^. (45) 

Due to this block structure, we can restrict the search 
over Z in the dual program, to Z that have the same 
structure, so we can take 

z = Zoe^i^-* ezj«, (46) 

where the Zi are operators in TCa ® T^-b ^ "Ha- The posi- 
tivity condition on Z in (|44|l . translates into a positivity 
requirements for each of the blocks in 1)46(1 . Using this 
structure we can write 

TtlFoZ] = Tt[Go{Zo + Zi + Z2)], (47) 

since Tr[Go^^Zf^] = Tr[GoZi], for i = 1,2 and X = 
A, B. We defined Gq in equation H19|l as a linear function 
of p, so we can write Gq = A(p), where A is a linear map 
from operators on Ha ^Hb io operators on Ha ®Hb® 
Ha, whose action on an arbitrary operator Y on Ha®Hb 
is given by 

k{Y) = Y(g)lA/dA + P(Y(E)lA)P/dA 

-]l^®TrA[r] 1^/4, (48) 

where P is the swap operator defined by P\i) ^\k)^ \j) = 
|j) (8) \k) (g) We can now define an operator Z on 
Ha'^Hb given by Z = A*{Zo + Zi + Z2), where A* is the 
adjoint map of A and is defined as the map that satisfies 
Tr[A(A:)y] = Tr[A:A*(y)] for any hermitian operators 
X, Y. For our particular case, this map takes the form 

A*iV) = TTc[V]/dA + TTc[PVP]/dA 

-lA®TTAc[V]/dl. (49) 

Then we have 

Tr[pZ] = Tr[A(p)(Zo + Zi + Z2)] = Tr[FoZ]. (50) 

Let psep be any separable state. Then we know that 
there is a PPT symmetric extension of Psepi or equiva- 
lently, the primal problem H23|l is feasible. Then from 
(|13|l . and using the fact that c = 0, we have that 
Tr[FoZ] > for aU dual feasible Z and so from ^ 
we have 

Tr[psepZ] > 0, (51) 

for any Z obtained from a feasible dual solution Z . This 
means that any operator Z constructed in this way, sat- 
isfies one of the two properties required in (|41|l . and is 
therefore a candidate for an entanglement witness. 



Now consider the case in which the primal problem is 
not feasible for a given state p. This can only occur if 
this state is entangled. We can then use the arguments 
presented in appendix IbI to affirm that there must be a 
feasible dual solution Zew that satisfies Ti'IFqZew] < 0. 
Using ((50|l and ((51|l we can see that the corresponding 
hermitian operator Zew satisfies the two conditions 

Tt[pZew] < and Tr[psepZEw] > 0, (52) 

which means that Zew is an entanglement witness for 
the state p. 

Even though we have shown the calculation explicitly 
only for the second test of the hierarchy, similar reasoning 
can be applied to all tests to show that if the primal 
problem is infeasible, there is a dual feasible solution that 
can be used to construct an entanglement witness for the 
state p. The EWs obtained for each of the tests have very 
well-defined and interesting algebraic properties, that can 
also be used to interpret each step in the hierarchy as a 
search for EWs of a particular form. 



A. Algebraic properties of the entanglement 
witnesses 

For any EW there is an associated bihermitian form 
given by ((43|l . We have shown that the requirement 
that an entanglement witness W is positive on all sep- 
arable states, is equivalent to requiring the associated 
form E{x, y) to be positive. 

Let us consider an EW obtained form the first test 
in the hierarchy, which corresponds to the usual PPT 
criterion. It is a well-known result that all states that 
fail this criterion, can be shown to be entangled by an 
entanglement witness of the form 

W^ = P + Q^-^, (53) 

where both P and Q are positive semidefinite operators. 
Entanglement witnesses that have this form are called 
decomposable. If we note by \tpp) the eigenvectors of P 
and by \<f>p) the eigenvectors of Q, we can write 

P = ^'^p\'^p){''Ppl 
p 

Q ^\Mp){(t>pV 
p 

where the eigenvalues Kp and Ap are nonnegative, since 
both P and Q are PSD. If we study the associated form 
Ew{x,y), we have 
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p p 

P ij P ij 



With l^p) = E.,>f,Kj) and = E.,<^r,Kj)- The 
last equality in (|54|l shows that Ew{x,y) can be written 
as a sitm of squared magnitudes (SOS), which proves its 
positivity. This property is an alternative description of 
decomposable entanglement witnesses. 

Now imagine that we have a state p that is PPT entan- 
gled, whose entanglement is detected by the second test 
of the hierarchy (i.e., p does not have a PPT symmet- 
ric extension to two copies of party A). Then we know 
that the dual SDP will provide us with an entanglement 
witness Zew for this state. Let us concentrate on the 
properties of this Zew- First, it is clear that it cannot 
be decomposable, since decomposable EWs can only de- 
tect states that are not PPT. By setting p^ep — \xy){xy\ 
in <|FT|l . we have that 

TY[\xy){xy\ZEw] = {xy\ZEw\xy) 

= Ez,Jx,y)>0. (55) 

According to equation H5U|) . we have 
TT[\xy){xy\ZEw] ^ TiiA{\xy){xy\){Zo + Zi + Z2)]. (56) 

The operator A maps a state p in Ha 'Si'Hb into an op- 
erator in Ha ®Ti.B 'Si'Ha that is invariant under swaps of 
the two copies of A and yields the original state p when 
one of the copies of A is traced out, but is in general not 
positive semidefinite. Now consider the state \xyx){xyx\. 
This state is invariant under swaps of copies of system 
A and also satisfies TTc[\xyx){xyx\] = \xy){xy\. Then 
we know that there must exist some coefficients aj such 



Since Zo^Zf"^ and are positive by construction, 

equation l|HT|l gives an exphcit sum of squares (SOS) de- 
composition of the right-hand side of (|60|l . We can con- 
clude then that even though the form E^^^^ (x, y) is not 
a SOS, it becomes a SOS when multiplied by the strictly 
positive SOS form {x\x) — \xi\'^. This property holds 
for any EW obtained from the second test. 

This result generalizes to all steps of the hierarchy: the 
bihermitian form associated with an EW obtained from 
the (fc-l-1)*'' test of the hierarchy, can be written as a SOS 



I 

that 

\xyx){xyx\ = k{\xy){xy\) + ^ajG j , (57) 

J 

since the Gj form a basis of the space of matrices M 
satisfying the swapping symmetry and Tr(7[Af] = 0. Ac- 
cording to (|44|l we have Tr[GjZi] = 0, and hence we can 
rewrite (|56|l as 

"TY[\xy){xy\ZEw] = "Tv[\xyx){xyx\{Zo + Zi + Z2)]. (58) 

Combining (|55|l and H58|) . we have 

{xy\ZEw\xy) = {xyx\{Zo + Zi + Z2)\xyx). (59) 

Since we are working with normalized states, we know 
that (a; I a;) = 1, so we can multiply the left-hand side of 
by this factor without changing the equality, obtain- 
ing 

= {xyx\{Zo + Zi + Z2)\xyx). (60) 
This equation is, in principle, only valid when the vari- 
ables Xi and yi correspond to a normalized state, i.e., 
when = 1 and Yli Ij/iP = 1- However, since both 

sides of (|6U|I are homogeneous functions of fourth degree 
on the Xi, and of second degree on the j/i, we can extend 
this equality to all values of the variables, and interpret 
H6(J|I as an equality between two forms that is satisfied 
everywhere. But we can now rewrite the right-hand side 
of ^ as 



(61) 

I 

when multiplied by the SOS form {x\x)'' — (X^il^jP)*^- 
We will say then that these EW are k-SOS. Then for 
example, an entanglement witness that is 0-SOS is de- 
composable, since its associated form can be written as 
a SOS (we will use SOS instead of 0-SOS for this par- 
ticular case). It is clear that if an EW is k-SOS, it is 
also l-SOS for all I > k. Note that for fc > 1, all k-SOS 
entanglement witnesses are indecomposable. 

As we discussed in Section V, searching for symmetric 
extensions with no PPT requirement generates another 



I 

{xyx\{Zo + Zi + Z2)\xyx) {xyx\Zi)\xyx) + {x*yx\Z'^'^\x*yx) + {xy*x\Z2''\xy*x). 



12 



complete family of separability criteria. For this family, 
equation H6U|) takes the form 

^ZB„-(2^'y)(2;|a;) = {xyx\Zo\xyx), (62) 

since now the LMI has only one block, corresponding to 
the positivity requirement on the extension. Since Zq is 
PSD, the right-hand side of still is a SOS, so we 
will still say that Zew is 1-SOS (or k-SOS if we replace 
{x\x) by (a; I a;)'' in (jHSJ)- The main difference between 
(|62|l and H61|l is in the type of terms that appear in the 
sum of squares decomposition. Note that (|62|) involves 
only squares of polynomials in the variables (xi, ijj) while 
the second and third terms in the right-hand side of H61|l 
correspond to squares of polynomials in the {x*,yj,Xk) 
and (xi, y*) variables respectively. This situation extends 
to all the steps of the hierarchy. The SOS decomposi- 
tion generated by the PPT family involves the squared 
magnitudes of all possible polynomials in the variables 
{xi,yj) and their conjugates that are compatible with 
the symmetry requirements, while the SOS decomposi- 
tion obtained from the non-PPT family involves squared 
magnitudes of polynomials involving only the variables 
{xi,yj). The completeness theorem tells us that for any 
EW obtained from any of the two families, we could make 
the associated bihermitian form into a SOS by multiply- 
ing by a certain power of the SOS form {J^i NiP)- How- 
ever, the value of the power needed in the non-PPT case 
will be in general higher than for the PPT case. 

We presented our family of separability criteria as the 
search for an extension of the original state p that sat- 
isfied certain symmetries and had positive partial trans- 
poses, and showed that this search could be put in the 
primal form of a SDP. From the discussion above, we 
can see that by looking at the dual SDP, we can inter- 
pret this hierarchy as a search over possible entangle- 
ment witnesses for the state p that are k-SOS for some 
k. Note that this interpretation also applies to the non- 
PPT family of criteria, the only change being the type 
of terms allowed in the SOS decomposition. The com- 
pleteness result proved in Section V tells us that the set 
of all entanglement witnesses that are k-SOS for some fc, 
is sufficient for proving entanglement of any state. This 
raises a very interesting question: are all EW k-SOS for 
some fc? We will now show that most EW are, and the 
only ones that may not be are those that are extremal in 
the sense that their associated hyperplane touches the set 
of separable states. These are the optimal EWs from [^. 

Let us first introduce some definitions regarding con- 
vex sets. A set K is said to be a convex cone if it is 
convex and closed under linear combinations with non- 
negative coefficients, i.e., ii x,y £ K and a, 6 > 0, 
then ax + hy G K. The dual cone of K is defined as 
K* — {z : {z,x) > 0,yx e K}, where (,) represents 
some inner product (note that K* will be different for 
different inner products). It is easy to show that K* is 
actually a closed convex cone even if K is not closed. An 
important property is that jl^ 

{Kr=d{K), (63) 



where c\{K) represents the closure of K. 

Let S be the set of all unnormalized separable states. 
5 is a closed convex cone. Its dual cone is S* = {Z : 
Tr[Zpsep] > 0,Vpsep G S}, which contains the set of all 
entanglement witnesses. If (5*)° notes the interior oi S* , 
we have {S*)" - {Z : Tr[Zp,ep] > 0,Vp,ep e S}. 

Theorem 2 Let W be an entanglement witness such that 
W G {S*)°. Then W is k-SOS for some k, i.e., 3k such 
that Ew{x,y){J2^ is aSOS. 

Proof: Let Ok = {Z : Ezix,y)iJ2^\x,\^)'' isaSOS}. 
This is just the set of entanglement witnesses that are k- 
SOS. Clearly, Ok C Ok+i and Ok C S* . Now we define 
the set 

oo 

O^lJOk. (64) 

O is a convex cone, although it may not be closed. We 
will now show that the dual of this cone is the set S. Let 
Psep G 5*. For any Z £ O, 3k such that Z G Ok. But 
Z G iS**, so TilZpsep] > 0, which means that psep G O* , 
so we have 

5 CO*. (65) 

Now, let p O* and assume p ^ S; then p is an en- 
tangled state. By the completeness of the hierarchy of 
separability tests, we know that there is a value of k for 
which Tr[Zp] < for some Z g Ok C O, and then we 
must have p ^ O* , which is a contradiction. Then 

O* c S. (66) 

^From H65|) and (|66|l . we have S — O*. Then we can use 
(|63|l to state that S* = cl(0), which means 

{S*y c O. (67) 

liW e (S*)" then by ^ there exists k such that W G 
Ok, and hence it is k-SOS. □ 

This theorem has a very nice geometric interpretation. 
It says that the sequence of convex cones Ok approxi- 
mates the convex cone of all entanglement witnesses S* 
from the inside, giving a complete characterization of its 
interior in terms the k-SOS property. On the other hand, 
the entanglement witnesses on the boundary of S* may 
not be k-SOS for any k. They satisfy Tr[Zpsep] — for 
some separable state psep , and correspond to the optimal 
entanglement witnesses discussed in 

As we briefly mentioned in Section II, the separabil- 
ity criteria based on searching for certain extensions of 
a state can easily be generalized to the study of multi- 
partite entanglement. The dual formulation of searching 
for an EW with certain algebraic properties also clearly 
applies to the multipartite case. It is known that mul- 
tipartite entanglement cannot be characterized in terms 
of bipartite entanglement alone However, our 

approach can be generalized to the multipartite case in 
order to construct another sequence of tests that is also 
complete. These results will be reported elsewhere |45j . 
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VII. EXAMPLES 

We present now some examples for which we ap- 
pUed our techniques to prove entanglement of certain 
PPT entangled states, and to construct the appropri- 
ate entanglement witnesses. For all these examples, 
the second test of the hierarchy (searching for PPT 
symmetric extensions to two copies of party A) was 
sufficient to show entanglement. We used MATLAB 
to code the corresponding SDP, and used the pack- 
age SEDUMI 46] to solve it. The code is available at 
http : //www. iqi . caltech . edu/publications .html. 

A. 3 ® 3 state. 

We consider the following state, described in 0, given 
by 

p„ = -|^+)(^+| + -a+ + ^Va+V, (68) 



with < a < 5, = ^ Elo N*)- = i(|01)(01| + 
|12)(12| + |20)(20|) and V the operator that swaps the 
two systems (note they are both the same space) . Notice 
that Pa is invariant under the simultaneous change of 
a — *■ 5 — a and interchange of the parties. The state is 
separable for 2 < a < 3 and not PPT for a > 4 and 
a < 1, which was proved in by using a positive map 
that is not completely positive due to Choi J20|. Our 
code solves the SDP for this state in about 5 seconds 
on a desktop computer. From this solution, numerical 
entanglement witnesses can be constructed for in the 
range 3 -I- e < a < 4 (and 1 < a < 2 — e) with e > 
10~^. A witness for a > 3 can be extracted from these 
by inspection: 



JEW 



2 (|00)(00| + |11)(11| + |22>(22|) + |02)(02| + |10>(10| + |21)(21| - 3|V+)(V^+|. (69) 



^From this entanglement witness, the Choi form 20] can be recovered. This observable is nonnegative on separable 
states: 

A{xy\ZEw\xy){x\x) = 3 \x2Xayl - xiX2ya\'^ + 3 \xixlyQ - X2XQy2\'^ + 3 \x2X2y2 - XQX2yaf + 3 \x2xly2 - 

+ \2xoXgy2 - 2xix*2yi + X2X2y2 - xoX2yo\'^ + l2a;oa;*yo - 2x2X2yi + X2xly2 - xixlyi\'^ 
+ \2xoXiy2 - X2Xoyl - xiX2yof + \2xoXQyo - 2xiXQyi + xix'lyo - X2xly2f' > 0. 



The expected value on the original state is Ti[ZEwPa] = 
^(3 — a), demonstrating entanglement for all a > 3. Ap- 
plying the non-PPT tests to this state fails to show en- 
tanglement for a < 3.84, even if we apply the 6*'* test, 
showing that this hierarchy can be considerably weaker 
than the PPT hierarchy. 

B. 4 (g) 4 state 

We consider next the 4 (g) 4 state given by ji^] : 

Pa = + \^2){ip2\ + a-a), a > 0, 

2 + a 

where 

V'l = i(l00) + lll) + \/2l22)), 

^2 = i(l01) + 110) + V2l33)), 

a = i(l02)(02l + 103)(031 + 112)(121 + 113)(131 -f 
8 

+ 120)(201 + 121)(211 + 130)(301 + 131>(311). 



Applying the PPT criterion yields provable entanglement 
only for those states with a < 2^/2 w 2.82843. It was sus- 
pected |47| that the state was actually entangled for all 
nonnegative values of a. Using our new criteria, we show 
that this is indeed the case, and provide an explicit en- 
tanglement witness and its decomposition. Again, only 
the second level of our hierarchy is needed. Using essen- 
tially the same approach as in the example above, from 
the dual solution of the semidefinite program we identify 
a particular witness: 



14 



W = (|22) - |00))((22| - (001) + (|22) - |11))((22| - (11|) + (|33) - |01))((33| - (01|) + (|33) - |10))((33| - (10|) 
+ |23)(23| + |32)(32| - |22)(22| - |33)(33|. 

This witness is nonnegative on all product states, as the following identity certifies: 



{xy\W\xy){x\x) = \xQxlyo + xixjyo - X2xly2 - xaXi^al^ + \xoxlyi + xixlyi - X2xly2 - X3X*Qy3f + 
\x2xly2 + X3xly2 - XQX^yQ - xix*2yi\^ + \x2X2y3 + x^x^y^ - xix^yo - XQxlyi\^ + 
\x1x3y2 - xoX2vl\^ + \xox3y2 - xiX2yl\^ + \x1x2yQ - xoX2vl\^ + |a;oa;32;o - xix?>yi? > 0. 



Applying the witness W to the state, we obtain 

Tr[VFpQ] = — ^ therefore certifying entan- 

glement for all values of a in the allowable range. 

We have applied the second test of the hierarchy to 
many bound entangled states found in the literature of 
dimensions up to 6 by 6 . In all cases, this test has been 
sufficient to demonstrate entanglement and construct nu- 
merical (and in some cases analytical) entanglement wit- 
nesses. However, we know from complexity arguments 
that there must be states that pass the second test in 
the hierarchy but are nonetheless entangled as we will 
discuss in the next section. 



VIII. PROPERTIES OF THE SETS OF 
ENTANGLED STATES WITH PPT SYMMETRIC 
EXTENSIONS 

For every k, let us consider the set of states that are not 
detected by the k*'^ test of the hierarchy, which can also 
be characterized as the states having PPT symmetric ex- 
tensions to k copies of party A. They generate a sequence 
of nested sets, each one containing the set of separable 
states. The completeness theorem tells us that this se- 
quence actually converges to the set of separable states. 
It is natural then to try to understand the properties of 
these particular sets. First of all, it is not difficult to see 
that these sets are all convex and compact. But there are 
other interesting questions we can ask about them. Are 
they nonempty? Is this sequence infinite or does it col- 
lapse to a finite number of steps for certain cases? What 
is the volume of the subset of entangled states contained 
in each set? Are these sets invariant under LOCC? In this 
section we will address some of these questions by explic- 
itly constructing states with PPT symmetric extensions 
and studying the implications of their existence. 

A. Constructing entangled states with PPT 
symmetric extensions 

We will now show explicitly how to construct an en- 
tangled state that passes the second test in the hierarchy, 
which means it has a PPT symmetric extension to two 
copies of system A. We will proceed by studying the 



properties of an EW obtained from the second test un- 
der a particular scaling and then use duality arguments 
to infer the existence of the required entangled state. The 
procedure is based on a scaling technique developed re- 
cently by Reznick 

Let p be a PPT entangled state that is detected by the 
second test of the hierarchy. Then we know that there is 
an entanglement witness Z that satisfies Tr[pZ] < and 

{xy\Z\xy){x\x) is a SOS, (70) 

which implies that TT[Zpsep] > for any separable state 
Psep- Since p is PPT, we know that the bihermitian form 
{xy\Z\xy) cannot be a SOS, otherwise Z can easily be 
shown to be decomposable and hence unable to detect a 
PPT entangled state, which would be a contradicition. 
Now, let us consider the following hermitian operator 

Z^ = {A-y (x)1b Z A-^ (E)1b, (71) 
with A — diag(l, 7, . . . , 7), 7 > 0. This operator satisfies 

{xy\Z^\xy) - {{A-'x)y\ Z\{A~'x)y) > 0, (72) 

since Z is positive on product states. Now consider the 

state 

p^ ^ j^{A(»lB)p{A^ <»1b), (73) 

with N a positive normalization constant. The state p-y 
is entangled for all 7 > 0, because 

Tr^Z^] = ^Tr[pZ] < 0. (74) 

This, together with H72|l . proves that Z~^ is actually an 
entanglement witness. Now, let us assume that Z-y is 
1-SOS for aU 7 > 0, that is 

{xy\Z^\xy){x\x) is a SOS. (75) 

By using 1)71(1 and introducing the variables x — A~^x, 
we obtain 

{S:y\Z\S:y){AS:\Ai) is a SOS, (76) 
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where the SOS structure is preserved due to the linearity 
of the transformation from x to x. Speciahzing for A in 
(f7?)|l we get 

{xy\Z\iy){{S:l+j'{S:l + ... + i^J) is a SOS. (77) 

Since the set of forms that can be written as a SOS is 
a closed set, the assertion in (f77)l must be valid also for 

7 = 0, and in that case we must have that 

{xy\Z\S;y)S:l is a SOS. (78) 

But this implies that the form {xy\Z\xy) is itself a SOS, 
which contradicts our assumption. Then we conclude 
that there must be a value 7* > such that for all 7 < 7*, 
the entanglement witness Zj has the property that 

{xy\Z-^\xy){x\x) is not a SOS. (79) 

Note that this argument is only based in the fact that 
there is a PPT entangled state detected by the second 
test of the hierarchy. Whenever this happens, we know 
that we can construct an EW that satisfies (|79|l . regard- 
less of the dimensionality of the subsystems. 

We can now very easily show that the existence of an 
EW satisfying H79|l implies the existence of an entangled 
state that is not detected by the second test. We just 
need a very simple lemma of convex analysis. 

Lemma 1 Let Ki and K2 be two closed convex cones 
such that Ki C K2, where C represents strict inclusion. 
Then the dual cones satisfy K2 C K*. 

Proof: Recall the definition of the dual of a cone K, 
K* = {z\{z,x) > 0,Va; e K). Let z € K^, then {z,x) > 
Oyx e K2. Since Ki C K2, {z,x) > 0,Vx G Ki, so 
z e K'^ and hence K2 C K'^. Now, let x G K2 such that 
X ^ Ki. Assume that K2 = K*. Let z G K^; then we 
also have z G K2 and so (z, x) > 0. But since this is true 
Vz G KI, this means that x G K^* = Ki, since Ki is 
closed, and this is a contradiction. Then we must have 

In our case we have the closed convex cones Ok = 
{Z : {xy\Z\xy){Y,^\x^\'^f is a S'OS'}, fc = 0, 1, .. ., that 
we have already defined for the proof of Theorem 2. In 
this section we have shown that if there is an entangled 
state that is detected by an EW in Oi that is not in 
Oq, then Oi is strictly contained in O2. According to 
the lemma above this means that O2, which is the set 
of states that are not detected by the third test of the 
hierarchy, is strictly contained in 0\, the set of states 
that are not detected by the second test. Thus, that 
there has to be a state that is not detected by an EW in 
Oi, because it belongs to 0\ (equivalently, it passes the 
second test), but is detected by an EW in O2, (because 
it does not belong to O2) and hence is entangled. 

The discussion above shows that there has to be an 
entangled state that passes the second test, but it does 
not give an explicit construction. But since we know that 
Z^ ceases to be in Oi for some small value of 7, and since 



we have shown that it detects the entanglement of p-y for 
all 7 > 0, we could be tempted to say that then p-y has 
to be a state that passes the second test for small enough 
7. The problem is that even though Z-^ is an EW for p^, 
it is not clear that there is not another EW in Oi that 
detects the entanglement of p-y. However, even if we are 
not assured that it would be the case, we can still check 
whether actually passes the second test for small 7. 
We did this numerically using our code for the case of 
the Choi state H68I) with a — 3.0001, and we found that 
indeed there is a value 7* ~ 0.4901 of the parameter such 
that for all 7 < 7*, the state p^ is entangled but cannot 
be detected by the second test of the hierarchy. 



B. Properties of the hierarchy 

^From an algebraic point of view, the result of the 
previous subsection is related to the fact that the fixed 
multiplier approach to proving nonnegativity by finding 
a SOS decomposition does not behave well under linear 
transformations. A solution to this problem may be al- 
lowing the multiplier {x\x) to vary as well, although it 
appears that this approach to checking for entanglement 
cannot be stated as a SDP. 

^From a physical point of view, this result has a very 
interesting interpretation. The transformation repre- 
sented by Equation (|73l) corresponds to applying an el- 
ement of a POVM that acts locally on system A and 
leaves system B alone. Such a transformation can be 
implemented by local operations with some finite prob- 
abihty. We see then that by SLOCC (stochastic local 
operations and classical communication), we can trans- 
form a state that is detected by the second test into a 
state that is not. Moreover, since the matrix A in H73|) 
is invertible, the reverse transformation is also possible 
under SLOCC. Then we could start with the state p', 
whose entanglement is not detected by the second test, 
and by LOCC operations obtain, with some probability, 
a state p that is detected by the second test. This shows 
clearly that, unlike the PPT class of states, the classes 
of states derived from the second and higher tests of the 
hierarchy are not invariant under LOCC. 

This scaling behavior of both states and entanglement 
witnesses is very general and has very important con- 
sequences on the hierarchy of tests. First, note that if 
we assume that Z^ is k-SOS (for any fixed k) for all 
7 > 0, this will be equivalent to replacing {x\x) by {x\x)'^ 
in H75|l . We can then follow the exact same steps dis- 
cussed after (|75|) and arrive to the same contradiction 
(i.e., that {xy\Z\xy) is a SOS). Then for 7 small enough, 
{xy\Z^\xy) {x\x)^ must cease to be a SOS for any fixed 
value of k. By applying Lemma 1 again, we can conclude 
that if there exists a PPT-entangled state p in T-La^'Hb, 
then for any value of k, there must he an entangled state 
Pk that is not detected by the fc*'' test of the hierarchy. 
Note that this result depends only on the existence of at 
least one bound entangled state, and not on the dimen- 
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sions of the system. Since there are expUcit examples of 
bound entangled states in 3® 3 (like Equation H68|) ') and 
2 (g) 4 (see ^13] ), and we can (by embedding) use them to 
construct bound entangled states in M, N x M > 6, 
we can conclude that there are always entangled states in 
Ti-A ^Ti-B, d-A X ds > 6, that pass the first k tests in the 
hierarchy, for any fixed k. In other words, the hierarchy 
never collapses to a finite number of steps, even for fixed 
dimensions (except in the already known cases of 2 (g) 2 
and 2 (8> 3). 

A very interesting question has to do with what is ac- 
tually the volume of the set of entangled states that are 
detected by the k*^ test, but are not detected by the 
{k — test (or equivalently, the set of states that have 
PPT symmetric extensions to (fc— 1) copies of A, but not 
to k copies). Even though we cannot give any estimate 
on the value of this volume, we can assert that this vol- 
ume is finite, i.e., the set in question has nonzero measure 
when we consider the measure on the set of states intro- 
duced in This is true for all values of k for which 
this set is nonempty. The proof of this fact is a straight- 
forward translation of the result presented in Lemma 7 
of ||4^ . This lemma proves that if a convex set Ci strictly 
contains a compact convex set C2 that itself contains a 
nonempty ball, then Ci contains a nonempty ball that 
does not intersect C2. In our case we can take Ci to 
be the set of states with PPT symmetric extensions to 
{k — 1) but not to k copies of A, and C2 the ones with 
extensions to k copies. Since both Ci and C2 are con- 
vex and compact, and C2 contains the set of separable 
states that contains a nonempty ball, the lemma proves 
that there is a nonempty ball of states that have PPT 
symmetric extensions to (k — 1) copies of A, but not to 
k copies. 



IX. CONSTRUCTING BOUND ENTANGLED 
STATES FROM INDECOMPOSABLE 
ENTANGLEMENT WITNESSES 

In previous sections we have discussed how to use 
semidefinite programs to implement separability criteria, 
and in particular we showed how to exploit the duality 
of the SDP to generate an indecomposable entanglement 
witness from a bound entangled state. In this section we 
will show that we can also use a SDP to test whether a 
given entanglement witness is decomposable. If the EW 
is indecomposable, the dual program constructs a bound 
entangled state that is detected by the witness. The re- 
sults of this section were reported in 50] which we follow 
closely. 

As discussed above, a sufficient but not necessary con- 
dition for any hermitian matrix Z to be an entanglement 
witness is for it not to be positive but rather decompos- 
able as Z = P + Q^-*, where P > 0, Q > 0. Such entan- 
glement witnesses are obtained whenever a state fails the 
first test of the hierarchy, which is just the PPT criterion. 
These entanglement witnesses can only detect entangled 



states that have a non-positive partial transpose. As it 
was shown in (|54ll , the bihermitian forms associated with 
them can be written as a SOS. 

If we know only the matrix elements of Z it may not 
be clear how to determine whether Z is decomposable or 
not. Consider the following semidefinite program in the 
dual form 

maximize — Tr [P + Q'^^]/d Ads 
subject to P > 0, Q>0 

H{P + QT^) = H{Z), (80) 

where H{Y) — Y — (Tr[y])l/(i^ds is a linear map that 
outputs the traceless part of Y. We can make (|Sn|) take 
the more familiar form H12I) if we introduce the matrix 
variable X, defined by A" — P ® Q {X will play the 
role of Z in (|12|l 'l. The final equality constraint may 
be enforced by a finite number {d\d^ — 1) of trace con- 
straints that define the matrices Fi and the coefficients 
Ci. Fq is then proportional to the identity. By adding a 
sufficiently large multiple of the identity to any matrix 
satisfying the trace constraints, it is always possible to 
construct anX = P©Q>0 also satisfying the con- 
straints. This means that the optimization is strictly 
feasible. Let be the optimum value of the objective 
function — Tr [P -|- Q'^-^J/d^ds, and Popt, Qopt be the 
values of P and Q that achieve this optimum. If 

V > -Tr [Z/dAde] , (81) 

then for e = TT[Z/dAdB] -f ry > it is clear that we can 
write 

^ = (Popt+elAB) + QrA' (82) 

which shows that Z is decomposable. 

We have stated the semidefinite program in its dual 
form. The primal form is worth considering since in the 
case where Z is nondecomposable, it constructs bound 
entangled states that are detected by Z. Using the for- 
mulae Hll|) , H12|l and H80|) , the primal form may be shown 
to be 

minimize Tr [Zp] — Tr [Z/dAds] 
subject to p > 0, p-^-* > 

Tr[p] = l, (83) 

where now the variables are the components of the state 
p in some basis. The completely mixed state is strictly 
positive and is a feasible solution of (|83|l . Thus, because 
of the strict feasibility of both the primal and dual prob- 
lems the optima of these two programs are equal 
and there are matrices popt , Popt , Qopt achieving the op- 
timum. By complementary slackness the range of Popt is 
orthogonal to the range of popt and the range of Qopt is 
orthogonal to the range of p^^^ . 

Suppose now that the optimum rj satisfies rj < 
—Tr[Z/dAdB]- This, together with rj — Tr[Zpopt] ^ 
Ti[Z/dAdB], means that 

Tr [Zpopt] < 0. (84) 
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Since we know that Z is an EW, then H84() means that 
the state popt is entangled. Furthermore, since this state 
is a feasible solution of (|S5|) . it must satisfy p^^^ > 0, so 
this state is bound entangled. For any P > 0, Q > we 
have 

Tr \[P + g^-)popt] = Tr \Pp,^,\ + Tr \Qpl;,\ > 0, (85) 

so Z cannot be decomposable, since it satisfies (|84ll . 

We will now show that popt is the so-called edge PPT 
entangled state. Since popt is a PPT entangled state, we 
can write popt = — p)Psep+pS, where psep is separable, 
(5 is a so-called edge PPT entangled state and p > 
is the minimum value for which such a decomposition 
is possible ,25j. An edge PPT entangled state 6 has the 
property that for any pure product state |x, y) and e > 0, 
d — ejx, y) (x, j/| is either not positive or not PPT. Since 
Tr \Zpscp\ > (because Z is an EW), if p < 1 then 

Tr[Zpopt] >Tr[Z(5]. (86) 

But (|86|l contradicts the optimality of popt i unless Popt is 
itself an edge PPT entangled state. 

This SDP finds the canonical decomposition of an 
indecomposable EW discussed in Defining e — 

— Tr[Zpopt] > 0, we have 

Z = Popt + Ql^t-^Ti-AB, (87) 

and as a result of the original dual form of the optimiza- 
tion, e is the smallest value for which such an expres- 
sion holds with P > 0, Q > 0. The range properties of 
^opt , Qopt and popt mean that this is the canonical form 
for Z introduced by Lewenstein et aL[25j|. 

X. CHARACTERIZATION OF POSITIVE MAPS 

It has been known for quite some time that there is a 
close relationship between entanglement witnesses, pos- 
itive bihermitian forms and positive maps lia. |2(]| |. In 
particular, this relationship was exploited in '11] to give 
a complete characterization of the separability problem 
in terms of positive maps. We will now show how to 
translate the properties of the entanglement witnesses 
generated by our hierarchy of separability tests into a 
characterization of the set of strictly positive maps. 

Let us denote by Aa and Ab the set of linear oper- 
ators acting on Ha and Hb respectively. We will call 
C{Aa,Ab), the set of hncar maps from A a to Ab- We 
say that a map A e C{Aa,Ab) is positive, if for any 
operator L e Aa, L > 0, then A(L) > 0. A completely 
positive (CP) map, is a map A such that the induced map 

An = A(E)K ■■ Aa<E) Mn ^ Ab 'E)Mn, (88) 

is positive for all n, with Mn being the space of opera- 
tors in a Hilbert space of dimension n and ]1„ the identity 



map in that space. CP maps have very important appli- 
cations in characterizing the set of physically meaningful 
evolutions of a quantum state. 

It is clear that any CP map is also a positive map. 
However, there are positive maps that are not CP. This 
has very important consequences on the study of entan- 
glement of quantum states. In particular, there is a one to 
one correspondence 11] between entanglement witnesses 
and positive non-CP maps. Since the hierarchy of separa- 
bility tests offers a characterization of the interior of the 
set of entanglement witnesses, it is not difficult to trans- 
late this characterization to the set of positive non-CP 
maps. To do this, we use the fact that for any linear op- 
erator L e Aa(!^Ab, we can define a map A £ £{Aa,Ab) 
by 

{k\A{\^){J\)\l) = {^\^{k\L\J)^\l). (89) 

Conversely, equation (|89|l can be used to uniquely con- 
struct the operator L from the map A. Equivalently, we 
can write [l^ 

A(p) = Tr^ [L{p^ ® 1b)] , (90) 

were p is an operator in Aa- Note that the same oper- 
ator L G Aa ® Ab can be used to define two different 
maps in C{Aa, Ab) and in C{Ab,Aa)- It was shown 
that this relationship gives in fact a one to one corre- 
spondence between entanglement witnesses, i.e., hermi- 
tian operators that are positive on separable states but 
have a negative eigenvalue, and positive non-CP maps. 
By using H9U|I it is not dificult to see that the interior 
of the set of entanglement witnesses, which correspond 
to those Z that satisfy Ti[Zpsep] > for any separable 
state Psep, is mapped onto the set of positive maps that 
map any nonzero positive semidefinite operator into a 
positive definite operator. Our characterization of entan- 
glement witnesses will translate into a characterization 
of this subset of positive maps. The maps that are left 
out are those that send at least one PSD operator into a 
another PSD operator that is not positive definite. 

In Section VI we showed that any Z in the interior of 
the cone of all entanglement witnesses is k-SOS for some 
k. Since they correspond to strictly positive maps (the 
ones that map any nonzero PSD operator into a positive 
definite operator), we can characterize these maps by as- 
sociating a bihermitian form directly to the map, using 
equation H89|) . Then we can state that a map is strictly 
positive only if the form 

EA{x,y) = {y\A{\x*){x*\)\y) 

ijkl 

is k-SOS for some value of k. 

We can also give an interpretation of this characteri- 
zation in a language that only involves statements about 
maps. To do this we need to analyze in more detail some 
of the properties of the EW generated by the SDP. Let us 
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consider the family of separability criteria that searches 
for symmetric extensions of a certain state, but does not 
require positive partial transposes. It is not difficult to 
see that the entanglement witnesses generated by the sec- 
ond test will satisfy 

{xyx\{ZEw ® tA)\xyx) = {xyx\Za\xyx) , (92) 

for all states \x) and \y), with Zq > 0. This is the ana- 
logue to equation (|60() . It is not difficult to show that 
this equality implies that the operators Zew ® and 
Zq actually coincide when they are restricted to the sym- 
metric subspace of the copies of system A. Furthermore, 
this is true for any number of copies of system A. If we 
denote by tt^ the projector onto the symmetric subspace 
of nf" (which we will denote by TC^/), we have 

7rfc(8)lB(Z_EW<XilA»('=-i))^fc'8ilB = (7rfe(g)lB)Zo(7rA;(g)lB). 

k 1 ^^^^ 

Since Zq is PSD on the space TLa ® Ti-B ® Ti!^^ ^\ its 
restriction to the tensor product of l-Cj^ and H. b remains 
PSD, which is the right-hand side of H93|l. The complete- 
ness theorem of Section V then tells us that if Zew is a 
strictly positive entanglement witness, then there must 
exist a finite k for which equation 193() is true. 

We can now use the isomorphism defined by (|89|l to re- 
state (|93ll in terms of properties of maps. First we use the 
fact that this isomorphism gives a one to one correspon- 
dence between PSD operators L £ Aa'^Ab and CP maps 
A G C{Ab,Aa)- Let A : Ab — > Aa be the positive non- 
CP map associated with Zewi and let A^ : Aa -^-h^jc 

be defined by Afe (p) = 7rfe(/9® l^»(fc-i))7rfc. Equation (jHOI) 
can be used to check that the map associated with the 
operator tt^ (g) 1b{Zew ^ Ia^c-i) Is is given by 

(Afe o A) : -> A^vk . (94) 

But since the right-hand side of H93(l is PSD, this map 
has to be completely positive. 

On the other hand, if A is not a positive map, then 
the map (A^ o A) cannot be completely positive for any 
k. This is true because the map A^ always maps a non- 
PSD matrix into a non-PSD matrix, as we can easily 
show. Let \i) be an eigenvector of a non-PSD operator a 
in Ha, with negative eigenvalue. Then (i|cr|i) < 0. For 
any k the vector belongs to the symmetric subspace 
H)^'' and satisfies TTfeji)®*^ = ji)®*^. Then we have 

((*!«'=) Afc(a) - ((z|®^-)7rfc(a®l|'=-i)^fe 

= ((i|®'=)(a®lf'=-i)(|z)^'=) 

{i\a\i) < 0, (95) 

and so Ak{a) cannot be PSD. Thus, we have the following 
result: 

Theorem 3 // the map A : Ab ^ Aa is strictly posi- 
tive, then there is a finite k such that the map (A^ o A) : 
Ab ■^H\'' ^s completely positive. If for some k the 
map (Afe o A) is completely positive, then A is a positive 
map. 



k 


a 


1 


0.4 


2 


0.58769 


3 


0.68556 


4 


0.72727 


5 


0.77663 


6 


0.80766 


7 


0.823529 


8 


0.846137 



TABLE I: Number of extensions and optimal value of a. 



Since this characterization of positive maps does not 
require solving a SDP, because we only need to check 
positivity of a matrix, it is interesting to study how ef- 
ficient this approach is in actually proving positivity of 
a map. To answer this question we consider the follow- 
ing example based on the case of the 3 (g) 3 state con- 
sidered in Section FVII Al Let the map Aq be defined as 
Ao(p) — —Tr[p]l3, where I3 stands for the identity map 
in 7^3. The map Aq lies in the interior of the cone of 
positive maps. Consider now a convex combination of 
Aq and the positive map A^^^, induced by the witness 
in equation (|69|l . i.e., 

A„ = (1 -a)Ao + aA^^^,, < a < 1. 

The map A^^^^ is in the boundary of the cone of positive 
maps. We have normalized the maps so that Aa(l3) = 
I3. Since for a = we have Aq, = Aq and for a = 1 
we have A^ = A^^^, the maps Aq, are contained in a 
line segment with endpoints near the center and in the 
boundary of the cone of positive maps, respectively. This 
implies that Aq is a strictly positive map for a < 1. 

A natural question in this case is to determine the 
ranges of a for which we can effectively recognize pos- 
itivity by applying the result of Theorem El For this, 
as explained, we have to form the tensor product of the 
given map with k — 1 copies of the identity, project on 
the symmetric subspace, and check whether the resulting 
matrix is positive semidefinite. The computation of the 
optimal a can be done in this case by solving a simple 
generalized eigenvalue problem. 

We have solved this numerically, for values of fc up to 8 
(this involves matrices of size 3 • i.e., 135 x 135). 

The obtained extreme values are shown in Tabled where 
k is the number of extensions. The results are consistent 
with the expected behavior limfe_>oo ctfe = 1. 

Notice that the convergence appears to be relatively 
slow, of order 1/fc; in contrast, the SDP tests presented 
earlier based on the PPT hierarchy can get all the way 
to the boundary a = 1 in just one step. 

It is interesting to note that Jamiolkowski also studied 
the problem of checking positivity of maps [IJ. His ap- 
proach was related to ours in the sense that he showed 
that checking positivity of a given map was equivalent 
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to the nonnegativity of a certain associated real polyno- 
mial. He then applied a general technique for checking 
positivity of polynomials. As discussed in the introduc- 
tion there are several such algebraic methods and they 
all scale badly with the problem size. In our case the 
specific problem of checking positivity of a linear map 
between matrix algebras has been reduced to a series of 
tests of matrix positivity, but none of them succeeds uni- 
formly for all maps. However, it is still the case that for 
many instances, the positivity of a given linear map can 
be determined and certified efficiently. 



XI. CONCLUSIONS AND DISCUSSION 

In this paper we have discussed a new family of sepa- 
rability criteria for bipartite mixed states. Each criteria 
consists in searching for an extension of a given state in 
a bigger space formed by adding a number of copies of 
one of the subsystems, and requiring this extension to be 
symmetric under exchanges of the copies and to remain 
positive under any partial transpose. A failure to find 
such an extension proves entanglement of the state, since 
it can be explicitly shown that separable states have the 
required extensions. If an extension is found, the test 
is inconclusive. This family of tests can be arranged in 
a hierarchical structure, with each test being at least as 
powerful as all the previous ones, and with the first test 
corresponding to the well-known Peres-Horodecki PPT 
criterion. 

This hierarchy of tests has two main properties that 
make it useful and appealing. First, the hierarchy is com- 
plete: any entangled state will fail one of the tests at some 
finite point in the sequence. Second, each test can be cast 
as a semidefinite program, which can be efficiently solved. 
Furthermore, by exploiting the dual structure of semidef- 
inite programs, whenever a state is proven to be entan- 
gled by failing one of the tests, an entanglement witness 
for that state can be explicitly constructed. This duality 
can also help us to interpret the hierarchy as trying to 
prove entanglement of a state by searching for entangle- 
ment witnesses with a particular algebraic property that 
states that the bihermitian form associated with the en- 
tanglement witness can be written as a sum of squares 
when multiplied by a fixed sum of squares to a certain 
power. The completeness of the hierarchy can then be 
used to show that this algebraic property characterizes 
all the elements in the interior of the cone of entangle- 
ment witnesses. 

We analyzed the computational resources needed to 
implement these tests. We found that for a fixed test in 
the hierarchy, they scale polynomially in the dimensions 
of the state. When we keep the size of the state fixed, 
the resources also scale polynomially with the number of 
copies added, or equivalently, with the order of the test 
in the hierarchy. This behavior is very interesting in light 
of recent results on the worst case complexity of the sep- 
arability problem. It has been shown that checking sepa- 



rability of a state is an NP-hard problem when we study 
the scaling with respect to the dimensions of both par- 
ties, so computational resources to solve it cannot scale 
polynomially in this general case. In our family of tests 
this non-polynomial behavior is reflected in how high up 
the hierarchy we need to go to detect all entangled states. 
Even though each test is efficiently implementable, there 
are states for which we need to go arbitrarily high in the 
hierarchy to show that they are entangled. 

The dual formulation of the hierarchy can also be un- 
derstood as the construction of a sequence of cones, each 
one containing the previous ones, that approximate the 
dual of the cone of separable states (which contains the 
entanglement witnesses) from the inside, giving a com- 
plete characterization of its interior. 

We can also interpret the primal formulation as the 
construction of a sequence of nested cones that approx- 
imate the cone of separable states from the outside. It 
is worth comparing this point of view with the results 
in 51], where a semidefinite program was used to approx- 
imate the cone of separable states from the inside. This 
result however, only applies when one of the subsystems 
has dimension 2, and gives a complete characterization of 
separability only in this particular case, while our hier- 
archy works for arbitrary dimensions of the subsystems. 

The hierarchy of tests allows us to divide the set of en- 
tangled states into different classes, according to whether 
they have PPT symmetric extensions to k copies of one 
of the parties or not. This generates a nested sequence of 
subsets of entangled states. This sequence can be shown 
to be infinite for all dimensions of the subsystems, ex- 
cept for 2 ® 2 and 2 (g) 3 where it is well-known that the 
PPT criterion is enough to characterize entanglement (in 
these two special cases, the hierarchy collapses to the 
first step). Furthermore, if the set of states with PPT 
symmetric extensions to (fc — 1) copies of A but not to k 
copies is nonempty, then it can be shown to have nonzero 
measure. These classes of states however, are not closed 
under LOCC operations, since we can transform a state 
that is not detected by the second test into a state that 
is, with finite probability and by applying only local op- 
erations. 



XII. ACKNOWLEDGEMENTS 

It is a pleasure to acknowledge stimulating conversa- 
tions with Hideo Mabuchi, John Doyle, John Preskill and 
Ben Schumacher. Thanks to Patrick Hayden for suggest- 
ing an improvement in an earlier proof of Theorem 1, 
and to Nicolas Gisin for providing us with the PPT en- 
tangled state in Section VII. B. ACD gratefully acknowl- 
edges conversations with Barbara Terhal. FMS thanks 
Oscar Bruno for many clarifying discussions regarding 
the completeness theorem. PAP acknowledges interest- 
ing conversations with Bruce Reznick. This work was 
supported by the National Science Foundation as part 
of the Institute for Quantum Information under grant 



20 



EIA-0086083, the Caltech MURI Center for Quantum 
Networks (DAAD19-00-1-0374) and the Caltech MURI 
Center for Uncertainty Management for Complex Sys- 
tems. 



APPENDIX A: IMPROVED SDP FOR 
IMPLEMENTING THE TESTS 

We will now introduce a slight modification of the SDP 
given in H23|l . that has the advantage of performing better 
numerically. With F given by (|22ll . let us consider the 
following SDP 

minimize t 

subject to Haba + F{x.) > 0, (Al) 

where Iaba is the identity matrix on the space Ha ^ 
I-Lb ® 'Ha- It is clear that we can always choose t such 
that the LMI on the second line of ljAl|l is satisfied. If 
the minimum of t is negative or zero, then there exists 
a value of x such that F(x) > 0, which is equivalent 
to say that (|23|) is feasible. On the other hand, if the 
minimum of t is strictly positive, then we know that -F(x) 
cannot be PSD. Thus we see that feasibility of (|23|l is 
equivalent to whether the minimum of ljAl|l is strictly 
positive or not. So we can use IjAip to detect entangled 
states. This approach has the property that the SDP 
(|A1(I is always feasible. This property makes the SDP 
solver to behave better numerically (because it uses an 
interior point algorithm). This is in fact the SDP that 
our code is solving when applying the tests to a given 
quantum state. 

APPENDIX B: STRONG DUALITY AND SDP 
INFEASIBILITY 

We want to obtain infeasibility witnesses for the SDP 
Fo + ^x,i^, >0. 

i=l 

Clearly, if we can find a lU > such that 
Tr[F,IU] = 0, Tr[FoIU] < 0, 



then the SDP is necessarily infeasible, as follows by the 
argument given after equation H14I) . Under what condi- 
tions does such a W exist? As we mentioned earlier, we 
need some form of strong duality to hold. 

Consider the set S := range -|- K, where K is the 
PSD cone, T : M™ §" is the linear map defined by 
•^W := E" 1 and A + S = {y\y = a + b,aeA,be 
B}. Feasibility of the SDP is equivalent to Fq e S. The 
set S is obviously convex. Now, if S is also closed, then 
we can apply the separating hyperplane theorem, and 
conclude the existence of a as above. 

The difficulty, of course, is than in general the sum of 
two closed sets may not be closed. In particular, in SDP 
things can go wrong. For instance, for 



which is obviously infeasible, it is not hard to see that no 
witness > as above can exist. This can be traced 
back to the fact that S in this case is not closed. 

So, what conditions can be required to guarantee that 
S be closed? An often-used criterion is the so-called 
Slater condition which in our case is the following: 

If kevT* n ri AT* 7^ 0, then range -|- K is closed. 

Here, JF* is the adjoint map of !F, K* is the dual cone 
(equal to K, in this case), and ri denotes the relative 
interior of a set. 

In other words, to guarantee the existence of infeasi- 
bility witnesses of the form we described (for any pos- 
sible Fq), it is sufficient to show a Z > 0, that satisfies 
Tr[FiZ] = 0, for all i = 1, . . . , m. Notice that this almost 
looks like the certificate W we are after, except that Fq 
does not appear in the expression (otherwise, the condi- 
tion would be useless). In general, checking whether the 
Slater condition is satisfied in concrete problems is not 
too difficult. 

For our SDPs in H23|l and H44|) , it is immediate to show 
that the criterion is indeed satisfied, as all the matrices 
Fi are traceless, so we can just take Z — 1 > 0. 
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